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^ We present photometric and spectroscopic data for the second 

microlensing event seen toward the Small Magellanic Cloud (SMC), 
MACHO-98-SMC-1. The lens is a binary. We resolve the caustic crossing 
o and find that the source took 2At = 8.5 hours to transit the caustic. We 

O 1 measure the source temperature T e g = 8000 K both spectroscopically 

£5 . and from the color (V — I)o ~ 0.22. We find two acceptable binary-lens 

models. In the first, the source crosses the caustic at <p = 43°. 2 and 
the unmagnified source magnitude is I s = 22.15. The angle implies 
^ ■ that the lens crosses the source radius in time = At sin = 2.92 

■ hours. The magnitude (together with the temperature) implies that 
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the angular radius of the source is 6* = 0.089 /xa,s. Hence, the proper 
motion is fx — 9*/t* = 1.26 kms -1 kpc -1 . For the second solution, the 
corresponding parameters are <fi = 30°.6, I s = 21.81, i* = 2.15 hours, 
9* = 0.104 //as, fx = 9*/t* = 2.00 kms -1 kpc -1 . Both proper-motion 
estimates are slower than 99.5% of the proper motions expected for 
halo lenses. Both are consistent with an ordinary binary lens moving 
at ~ 75-120 km s -1 within the SMC itself. We conclude that the lens is 
most likely in the SMC proper. 

Subject headings: astrometry, gravitational lensing, dark matter, galaxies: 
individual (Small Magellanic Cloud, Large Magellanic Cloud) 
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1. Introduction 

When microlensing searches toward the Large Magellanic Cloud (LMC) were 
initiated seven years ago, it was believed that the microlensing events themselves 
could resolve the question of whether the Galactic dark halo is composed primarily 
of baryonic material in the form of Massive Compact Halo Objects (MACHOs). The 
major potential difficulty was thought to be the problem of distinguishing lensing 
events from variable stars. 

The situation today differs markedly from these expectations. On the one hand, 
more than a dozen candidate events have been detected toward the LMC (Aubourg 
et al. 1993; Alcock et al. 1997a), and while a few may be intrinsic variables, the great 
majority are clearly due to genuine microlensing events. On the other hand, the 
interpretation of these events is far from clear. The optical depth is r ~ 2.5 x 10 -7 
(Alcock et al. 1997a), about half what would be expected from a full dark halo of 
MACHOs, and far more than could be due to known populations of stars along the 
line of sight in the Galactic disk, Galactic spheroid, and LMC disk. By contrast, 
from the measured Einstein crossing times £e ~ 45 days and the well-constrained 
kinematics of the halo, the mass of these objects would appear to be 0.4 M Q . These 
could not be made of hydrogen or the population would be discovered easily. The 
paucity of obvious solutions has led to a large number of alternative explanations 
for the observed events, including stars within the LMC bar/disk (Sahu 1994), a 
disrupted dwarf galaxy in front of the LMC (Zhao 1998, Zaritsky & Lin 1997), a 
warped and flared Milky Way disk (Evans et al. 1998), non-standard halo kinematics, 
and a halo composed of white dwarfs. However, all of these solutions present 
considerable difficulties (Gould 1995, 1998a; Beaulieu & Sackett 1998; Bennett 1998; 
Fields, Freese, & Graff 1998; Gyuk & Gates 1998). 

More observations of the type already in progress seem unlikely to resolve the 
issue by themselves; new methods and strategies will be required. One approach is 
to search for microlensing toward the Small Magellanic Cloud (SMC). If the LMC 
events are due to halo lenses, these should produce events of similar duration and 
frequency, though not identical if the halo is significantly flattened (Sackett & Gould 
1993). MACHO (Alcock et al. 1997b), EROS (Palanque Delabrouille et al. 1998) and 
OGLE (Udalski et al. 1998) are carrying out such observations but, since the SMC 
has many fewer stars (and hence fewer expected events), this strategy is unlikely to 
resolve the matter by itself within the next few years. 
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Another approach is intensive photometric monitoring of ongoing events 
with the aim of measuring either "parallax" or "proper motion." These measure 
respectively v/x and v/(l — x) where v is the transverse speed of the lens relative 
to the observer-source line of sight, and x is the ratio of the lens to source distance, 
and can distinguish between halo and LMC lenses almost on a case-by-case basis 
(Boutreux & Gould 1996). Typical events are long enough that the Earth's motion 
would often cause detectable "parallax" deviations in the light curves if the lenses 
were in the halo, provided that all events were monitored intensively for their entire 
duration (Gould 1998b). GMAN (Alcock et al. 1997c) is currently following all 
LMC and SMC events with the aim of detecting this effect, although probably with 
insufficient intensity to do so in most cases. EROS and MACHO have recently jointly 
proposed much more intensive follow-up observations at the European Southern 
Observatory, but have not as yet been awarded telescope time. Proper motions 
can be measured if a lens caustic (region of formally infinite magnification) passes 
over the face of the source. The finite size of the source suppresses this infinite 
magnification according to its size. By modeling the light curve, one can therefore 
measure t*, the time it takes the lens to cross the source radius. The angular size of 
the source, 9* can be estimated from its color, flux and the Planck law. The proper 
motion is then given by, 

V=— I 1 ) 

Unfortunately, the majority of events are caused by point-mass lenses for which 
the caustic is a single point. The probability that this point will pass over the face 
of a typical source is ^ 10~ 2 . However, for binary lenses of near-equal mass ratio, 
the caustics are closed curves whose size is of order the Einstein radius. This means 
that when there is one crossing, there is usually a second - with forewarning. 



On 8 June 1998, MACHO ( [http: / / darkstar.astro. Washington. edu| ) issued an 



alert for such a second caustic crossing for the second event discovered toward 
the SMC, MACHO-98-SMC-1 (J2000 a =00:45:35.2, 5 = -72:52:34). This alert 
provided a rare opportunity to determine the lens location for at least one SMC 
event. 
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2. Observations 

The Probing Lensing Anomalies NETwork (PLANET) was formed to monitor 
ongoing microlensing events on a round-the-clock basis from observatories girding 
the southern hemisphere with a view to detecting anomalies in the light curves which 
might betray the presence of a planetary companion to the lensing object (Albrow et 
al. 1998). We have obtained substantial dedicated telescope time during the bulge 
season (southern Autumn and Winter) on the SAAO 1 m (South Africa, 20° 49', 
-32°), the CTIO-Yale 1 m (Chile, 289° 11', -30°), the Canopus 1 m (Tasmania, 
147° 32', -43°), and the Perth 0.6 m (Western Australia, 116° 8', -32°). In addition, 
we fortuitously obtained additional observations from the CTIO 0.9 m. All of our 
previous monitoring was of bulge events because they are the most common and 
are concentrated in a relatively compact observing season. However, in view of the 
importance of MACHO-98-SMC-1, we focused a major effort upon it. 

Nightly observations (weather permitting) were made from the SAAO 1 m, 
the CTIO-Yale 1 m, the CTIO 0.9 m, and the Canopus 1 m from 9 June to 17 
June. These consisted typically of 5 to 20 min exposures depending on the source 
brightness, mostly in Cousins /, but with some in Johnson V. Reductions were 
performed at the telescope as detailed in Albrow et al. (1998). On 15 June, MACHO 
predicted that the second crossing would occur on 19.2 ± 1.5 June UT on the basis 
of a binary-lens fit which they put on their web site. From PLANET observations 
alone (and confirmed with the addition of earlier MACHO observations posted on 
their Web site) we were able on 17 June to refine the prediction of the peak to 
be 18.0 June UT and issued our own caustic alert including light-curve data. For 
the next two days, we devoted all available telescope time to this event, and then 
continued sparser observations into the beginning of July. 

In addition, we obtained a spectrum of the source during the caustic crossing on 
18 June when the source was I ~ 17 using the 1.9 m SAAO telescope and grating 
spectrograph at a resolution of about 7 A. 

3. Analysis 

The proper motion can be determined using equation ([!]) from measurements of 
the source crossing time, £*, and the angular source radius, The crossing time 
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is extracted from the light curve, while the source radius can be measured by two 
independent methods, spectroscopic and photometric. 

3.1. Measurement of 

A static binary lens with a uniform surface-brightness finite source generates 
a light curve characterized by nine parameters. Three of these are t^, to, and 
M o(> 0), the Einstein radius crossing time, the time of closest approach, and the 
impact parameter of a hypothetical event in which the binary is replaced by a single 
point-lens of the same total mass placed at the geometrical midpoint of the binary. 
Two other parameters are I s and lb, the magnitudes of the lensed source and of 
any unresolved (unlensed) light superposed on the source star, respectively. Three 
parameters specific to the binary character of the lens are the mass ratio q = M^jMx 
of the secondary to the primary (0 < q < 1), the projected binary separation d in 
units of the (combined) Einstein radius r E , and the angle a (0 < a < 2ir) between 
the binary-separation vector (M2 to Mi) and the proper motion of the source relative 
to the binary center of mass (see Fig. [I]). The direction of the relative motion 
between lens and source is chosen so that the lens system is on the right hand side 
of the moving source. Finally, p* = 0*/0e is the source radius in units of the angular 
Einstein radius. 

Once the full solution is found, the crossing time is given by £* = p*^E- In order 
to understand the sources of uncertainty in this determination, however, it is useful 
to write the equation as, 

U = At sin 0, At = p*tE esc 0, (2) 

where is the angle at which the source crosses the caustic (see inset to Fig. |l|). 
The quantity (2 At) is the time taken by the source to transit the caustic. If the 
caustic crossing is well monitored (as it is in this case) then At is determined to 
high precision. We find below that At = 4.25 hours. Most of the uncertainty in 
therefore comes from the angle 0, which is determined from the global fit of the light 
curve. If the data permit two or more discrete geometries, then sin can take on 
substantially different values. 

Figure |2] shows data (nightly binned except at the peak) together with the 
light curves for our models I and II. Figure [I] shows the lens geometry, the resulting 
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caustic structure, and the trajectory of the source relative to this structure for model 

I, which has parameters = 108.4 days, to — 27.83 May 1998 UT, uq = 0.065, 
I s = 22.15, I b = 21.43, a = 348°.2, q = 0.29, d = 0.59, and p* = 0.00112. In 
particular, = p*£e = 2.92 hours. The source center passed over the caustic at UT 
18.120 June. The overall fit has x 2 = 394 for 164 degrees of freedom. The model II 
parameters are t E = 80.3 days, t = 10.68 June 1998 UT, u = 0.036, I s = 21.81, 
I b = 21.89, a = 351°.4, q = 0.95, d = 0.56, and p* = 0.00112, with t* = 2.15 hours, 
and x 2 — 442. Note that the values of t* are quite different in the two models 
despite the fact that values of At (4.26 hours and 4.23 hours) are very similar. The 
difference is due to the different angles at which the source crosses the caustic, 

= 43°.2 and 30°.6 respectively. The differences in t* and I s lead to a proper motion 
that is larger for model II by a factor of 1.59. The model II x 2 is ~ 48 higher than 
for model I. However, the fact that model I has a better fit does not rule out model 

II, which gives an acceptable fit to the data. Hence, the difference in proper motions 
between the two models (which is much larger than the statistical error of each 
solution) gives a good indication of the systematic uncertainty. It is possible that 
this uncertainty will be removed as we acquire more data, or by combining our data 
set with those of MACHO and EROS. In particular, we note that model I appears 
more consistent with the pre-first-caustic data that MACHO put on their web site. 
However, at present we must allow that our proper motion could be off by a factor 
~ 1.6. 

We have also fit for V s and Vb holding all seven other parameters fixed at the 
values found for the / band solutions. For model I, we find [V, (V — I)] = (22.45, 0.3), 
(21.75,0.3), and (21.30,0.3) for the source, background, and total light respectively. 
For model II, we find (22.10,0.3), (21.80,-0.1), and (21.20,0.1). Initially, the MACHO 
web page estimated the total baseline values to be V = 21.8, V — R — 0.1. After 
the appearance of this Letter as a preprint, MACHO (A. Becker 1998 private 
communication) informed us that they had revised their estimate to V = 21 .4. In 
addition, OGLE (A. Udalski 1998 private communication) confirmed this value, 
finding V = 21.41 ± 0.23. Since our fits contain no data from the baseline, these 
independent measurements of the baseline flux serve as an important check on the 
viability of the models. 

Our determinations may be compared with those of the EROS collaboration 
(Afonso et al. 1998) based on a single night of data covering the end of the caustic 
crossing. They found a lower limit At > 3 hours and measured the end of the caustic 
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crossing at UT 7:08 ±0:02 18 June 1998. In our model I, At = 4.26 hours, and the 
end of the caustic crossing is at UT 7:11 18 June 1998. 

3.2. Measurement of 9* 

We estimate the angular radius of the source, from the light curve. The model 
I fit to the overall light curve yields an unlensed source magnitude / = 22.15. We 
determine the color of the source, V — / = 0.31 ±0.02, from its measured color near 
the caustic, transforming our instrumental magnitudes to standard Johnson/ Cousins 
bands using a sample of stars in our data set that were independently calibrated 
by OGLE (Udalski et al. 1998). Note that while the color of the source at any 
particular time on the light curve could be affected by blending, this effect is 
extremely small near the caustic because the source is magnified by a factor ~ 100. 
We assume an extinction Ay = 0.22 ±0.1 (see below). This implies I = 22.02 
and (V — I) = 0.22. At fixed (V — I) color, the effective temperature depends 
only very weakly on metallicity (which we take to be [Fe/H]= —0.75). We use the 
Yale Isochrones (Green, Demarque, & King 1987) to estimate T e g ~ 7900 ± 300 K. 
Combining this with the dereddened flux, and again using the Yale Isochrones, we 
find 9* = 0.089 ^as. 

We now investigate the effects of various errors. First an extinction error 
of 5 Ay, alters Iq by — 8Aj ~ — 0.65 Ay (Stanek 1996), leading to an incorrect 
estimate of the radius by (59*/ 9*) ~ (In 10/5)0. 65Ay. This would be more 
than compensated, however, by the fact that the color would be in error by 
5E(V — I) ~ 0A5A V . This in turn would lead to an incorrect estimate of the 
temperature and so of the / band surface brightness, 5(2.5 log S) ~ ^y5E(V — I), 
where 7 = [(dkiBy / 'dT) / '(dlnBj / 'dT) - and By and Bj are the 

Planck intensities at 0.55 /zm and 0.80 fim, respectively. This causes an 
error in the radius of (59*/9*) ~ —(In 10/5)0. A^5Ay. Thus, the net error is 
(59*/ 9*) ~ (In 10/5) (0.6 - 0.^)6 A v 0.25<L4 y , for T ~ 8000 K. 

By interpolating the extinction map of Schlegel, Finkbeiner, & Davis (1988) 
across the SMC, we find a foreground extinction of Ay = 0.12. We adopt an internal 
extinction of 0.1 ± 0.1 mag, and obtain Ay = 0.22 ± 0.1. From the discussion above, 
the extinction uncertainty yields an uncertainty in the radius of ~ 3%. The 300 
K uncertainty in the temperature, which arises primarily from the uncertain color 
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calibration of the models (M. Pinsonneault 1998 private communication) generates 
a 5% uncertainty in the radius. Finally, the ~ 10% uncertainty in the unmagnified 
flux of the source produces another 5% uncertainty in the radius, for a final estimate 
of 0* = 0.089 ± 0.007 pas. 

From the spectrum, we find that the source temperature is consistent with the 
value T = 8000 K that was estimated from the color (V — I) ~ 0.22, i.e., roughly 
an A6 star. The flux (Vo = 22.24) is in reasonable agreement with what one would 
expect from a metal-poor star at 8000 K (My ~ 3.2) at the distance of the SMC 
(~ 60kpc.) Figure [3] shows the source spectrum. 



4. Results 

Combining the measurements of and t* we find 

Q 

H = — = 1.26 ±0.10 km s _1 kpc _1 , (3) 

for model I, where we have assumed that the error in t* is much smaller than in 
9*. As noted in the previous section however, model II yields \x = 2.00 kms -1 kpc. 
Hence, the difference between acceptable solutions is much larger than the error 
in each solution. Figure f| shows the distribution of proper motions expected 
from a standard p*(r) oc r~ 2 isothermal halo characterized by a rotation speed 
v c = 220kms _1 , and assuming an SMC distance of 60 kpc (Barnes, Moffett, & 
Gieren 1993). Since the proper motion of the SMC is difficult to measure and not 
yet tightly constrained (e.g. Kroupa & Bastian 1997), we show curves for 4 separate 
assumptions: three with an SMC transverse speed of 250kms _1 oriented at 0°, 90°, 
and 180° relative to the reflex motion of the Sun, and the fourth with no transverse 
speed. In all four cases we find that the fraction of events with proper motions 
smaller than that given by equation (Q) is less than 0.1% for model I. For model II 
the fraction is less than 0.5%. This would seem to imply that the lens is not in the 
halo. 

The main alternative explanation for the event is that the lens is in the SMC. 
In this case, its transverse speed relative to the source is v ~ 75kms _1 for model I 
and v ~ 120 km s -1 for model II. The measured one-dimensional dispersion of the 
SMC carbon stars is a ~ 21 ± 2 km s" 1 (Hatzidimitriou et al. 1997). Figure ^ shows 
the expected distribution of proper motions assuming that both the lens and the 
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source transverse velocity distributions are characterized by this speed. That is, the 
SMC material is assumed to be characterized by a single component (but see below). 
The fraction of all events with proper motions fi > 1.26 kms -1 kpc -1 is ~ 9%. The 
fraction with fi > 2.00 km s -1 kpc -1 is < 0.1%. 

Another possibility is that we are seeing material that has been tidally disrupted 
and is moving transversely at 75-120 km s -1 . There is substantial evidence for such 
rapidly moving material in radial velocity studies of the SMC. The extensive HI 
observations of McGee & Newton (1981) show that the SMC has four principal and 
overlapping structures at velocities Vh = 114, 133, 167 and 192 kms -1 . McGee & 
Newton (1981) state that of the 506 lines of sight examined, 97% show evidence for 
multiple peaks. Staveley-Smith et al. (1997) present much more detailed HI maps 
and argue that these data imply not multiple components but rather a complexly 
structured and perhaps disintegrating SMC. Zhao (1998) has modeled the SMC as 
comoving tidal debris and predicted that such a structure should generate about 
one microlensing event per year. Torres & Carranza (1987) summarize evidence for 
4 peaks at Vh = 105, 140, 170, and 190 kms -1 from their own HII data plus earlier 
results for HI, planetary nebulae, supergiants, emission regions, and Call. However, 
the picture that they present is much more of correlations between radial velocity 
and position on the sky. Of course, to produce the observed proper motion, there 
must be material at two different velocities along the same line of sight. We therefore 
believe that it is plausible, but by no means proven, that there is material within 
the SMC that could produce the proper motion indicated by model II. (As discussed 
above, the proper motion associated with model I can be explained from the internal 
dispersion of clumps within the SMC alone.) 

If the lens does lie in the SMC, then the Einstein radius is approximately 
r E ~ (4:GMD LS /c 2 y/ 2 , where D LS is the lens-source separation. The observed time 
scale is = te/v = 108 days or 82 days for the two solutions. Hence, the mass is 

for model I and a factor 1.4 times larger for model II. Various studies have found 
that the SMC is extended along the line of sight. Although the original estimates 
of 20 kpc or more by Mathewson, Ford, & Visvanathan (1986) are probably too 
high, the depth may be as high as 10 kpc (Welch et al. 1987; Martin, Maurice, & 
Lequeux 1989; Hatzidimitriou, Cannon, & Hawkins 1993). Equation ([|) implies that 
if the lens-source separation is comparable to this depth, then the total mass of the 
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lens is that of a typical low-mass binary. The projected separation of the binary is 
dr-E = bvtv ~ 2.7 AU or dr^ ~ 3.1 AU, which are also quite reasonable for a low-mass 
binary. The period would then be P ~ 6(.Dls/5 kpc) 1//2 years, indicating that the 
effect of internal binary motion during the event is small. 

In sum, MACHO-98-SMC-1 appears to be consistent with lensing by a stellar 
binary in the SMC proper, and inconsistent with lensing by an object in the Galactic 
halo. 
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Fig. 1. — Geometry of model I. The diamond shaped curve is the caustic, i.e., the 
region of formally infinite magnification. The thick solid line and two thinner parallel 
lines indicate the source trajectory and the finite size of the source. The tick marks are 
in units of the Einstein radius crossing time, t-^ = 108.4 days. The two components of 
the binary are shown by circles whose relative sizes are proportional to their masses. 
Their projected separation in units of the Einstein ring is labeled d. The insert shows 
more clearly the finite size of the source as well as the angle <fi at which the source 
intercepts the caustic. 



Fig. 2.- Light curve of the PLANET data for MACHO-98-SMC-1. Shown are 
the data from the SAAO 1 m [circles), the CTIO 0.9 m (squares), the CTIO-Yale 
1 m (triangles), and the Canopus 1 m (asterisks). The inset covers about 0.6 days, 
corresponding to less than one tick mark on the main figure. The data are binned 
by day (except near the caustic) for clarity. The light-curve fits were done without 
binning. The two curves show model I (bold) and model II (solid). The two models 
have very similar crossing times, At ~ 4.25 hours (which are determined from the 
caustic crossing), but different angles at which the lens crosses the caustic (see Fig |l]) , 
= 43°. 2 versus = 30°.6 (which are determined from the overall light curve). The 
proper motion is determined from the combination \x = 9*/t* = 6**/(At sin 0) where 
9* is the angular size of the source. 



Fig. 3.— PLANET spectrum of MACHO-SMC-98-1 taken near 18.0 June 1998 UT. 
The obvious Balmer line absorption indicates the source is an A-type star, consistent 
with the spectral type estimated from the V — I color of the source. The strength of 
weak Call H and K absorption suggests the effective temperature of the star is close 
to 8000K, although the Balmer line absorption strengths seem to indicate a slightly 
higher temperature. Note that the spectrograph slit was not well aligned with the 
parallactic angle, so the overall spectral energy distribution may be affected. 
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Fig. 4. — Distribution of expected proper motions from Galactic halo (right) and 
non-tidal SMC (left) lenses. The SMC model is based on a one-dimensional dispersion 
of 21kms _1 for both the lenses and sources. The normalizations are set arbitrarily 
to unity at the peaks. Four models are shown for the halo, three where the SMC is 
moving transversely at 250 km s -1 at 0° (solid), 90° (dashed), and 180° (bold dashed) 
relative to the reflex motion of the Sun, and one where the SMC has no transverse 
motion (bold). The boxes shows the best fit and 1 a errors for the proper motion of 
the lens for model I (left) and model II (right). The central values lies below 99.9% 
(99.5%) of the halo distribution and above 91% (99.9%) of the SMC distribution. The 
observed value is inconsistent with the halo hypothesis. The lens and source could 
possess normal SMC kinematics or one of them could be part of SMC tidal debris. 
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